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There have been several experimental investigations of entanglement enhancement of a two-mode 
squeezed vacuum. In particular, conditional preparation by photon subtraction has been shown to 
improve correlations achieved with this entangled resource. Here we analyse the role of Gaussian and 
non-Gaussian measurement for entanglement concentration acting on a pair of two-mode squeezed 
states. We find stringent requirements for achieving further entanglement enhancement by a joint 
measurement setup on the two resources. 



I. INTRODUCTION 

In recent years, quantum optics experiments have 
demonstrated the possibility of quantum enhanced com- 
munication by coherent manipulation of states of light. 
Important and illustrative quantum communication tasks 
such as quantum cryptography [Ij, quantum dense cod- 
ing [2J and quantum teleportation [3j can be implemented 
within the continuous variables (CV) framework, by 
manipulating the quadrature amplitudes of electromag- 
netic modes, specifically in Gaussian states [IHB]. This 
type of implementation is readily achieved in the labo- 
ratory, relying on linear optics for deterministic trans- 
formations, and on the nonlinear process of parametric 
down- conversion (PDC) for generation of non-classical 
and entangled squeezed Gaussian states. A challenge 
for this approach is that the fidelity of CV quantum 
communication depends on the quality of entanglement 
and unit fidelities are only achieved by employing in- 
finitely squeezed states which are impossible to realise 
in practice. In addition, quantum transmission channels 
degrade shared entanglement, decreasing the fidelity of 
communication protocols. 

In order to improve the performance of CV quantum 
communication protocols, methods have been developed 
to enhance entanglement initially shared in the form of 
bipartite Gaussian states. The term entanglement distil- 
lation ^ serves as a label for protocols consisting of 
local operations and classical communication that con- 
vert a number of identical two- mode entangled states into 
a smaller number of more entangled two-mode states. 
More specifically, when the initial states are pure, the 
protocols are termed entanglement concentration. En- 
tanglement distillation cannot be achieved if both the 
initial states and the local operations are Gaussian [g- 

In practice, a non-Gaussian operation may be achieved 
by measurement with single-photon detectors [12j. Pho- 
ton subtraction and photon addition on the modes of 
an entangled two-mode Gaussian state have both been 
shown to increase entanglement. These operations have 
been extensively discussed in recent experimental [131 [14] 
and theoretical [IF] reviews. In the context of entangle- 



ment concentration, both multiple applications [16j and 
the coherent superposition [17J of annihilation and cre- 
ation operators - well-approximated by photon subtrac- 
tion and addition - have been considered. It has been 
shown that photon subtraction can improve the fidelity of 
teleportation with CV, while simple photon addition can- 
not [18-20J. Photon-subtracted states have been shown 
to be more tolerant to loss than Gaussian states [21). 
The enhancement of photon subtraction on the nonlocal 
behaivour of two-mode squeezed states has also received 
attention [22','23]. 

Experimentally, an enhancement of entanglement in an 
initial two-mode Gaussian state has been demonstrated 
by various methods: photon subtraction in either or both 
modes [24J, non-local photon subtraction [25j, and weak 
Gaussian measurement after the state has gone through 
a non-Gaussian noise channel [26j. 

In this work we compare different measurement oper- 
ations acting jointly on a pair of two- mode Gaussian en- 
tangled states with the aim of performing entanglement 
concentration. In particular, we explore the advantage of 
pairwise interaction of resource states in the first step of 
the iterative protocol introduced in [271 EH], employing 
photon subtraction as the entanglement-enhancing non- 
Gaussian operation. The approach consists of combining 
two photon subtracted states on beam splitters and per- 
forming a Gaussian measurement on two of the output 
modes. It represents an example of two-source entangle- 
ment distillation feasible within the current experimental 
state of the art. 

We find analytical expressions for the result of inter- 
fering different photon subtracted states. We construct 
the best Gaussian measurements required to combine the 
entangled states and identify the arrangement for achiev- 
ing the optimal trade-off between success probability and 
enhancement of the entanglement. Finally, we analyse 
the resilience to loss of the approach. Our results show 
that, if limited to the first step of the iterative distillation 
protocol, combining photon-subtracted states has limited 
advantages with respect to simple photon subtraction. 
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II. ENTANGLEMENT CONCENTRATION BY 
PHOTON SUBTRACTION ON A TWO-MODE 
SQUEEZED VACUUM 

A two-mode squeezed vacuum (TMSV) is commonly 
produced by spontaneous non-degenerate PDC, a process 
occurring when a noncentrosymmetric crystal is pumped 
by an intense laser field. A photon can be spontaneously 
annihilated in the pump field and one photon created in 
each of the two modes a and 6, resulting in the state 



|^rM5v) = v/l-A2^ A" In) Jn)j 



(1) 



Eq. |3j the measurement operator corresponding to the 
subtraction of k-photons from mode a is Mk,a 



Mk,a =a, {k\ B{0)aa. |0),^ = 

x^r t' 



(4) 



The unnormalized state created by subtracting one pho- 
ton from a TMSV is found to be 



|^i,o) = v/l-A2rA^(Atr\/^ |n)„ |n + 1), (5) 



where \n) ^ is a number state with n photons in mode c 
and A = tanh(r) is determined by the parameter r that 
is proportional to both the nonlinear coefficient of the 
PDC crystal and the intensity of the pump laser beam. 
A TMSV is the typical entangled resource for CV quan- 
tum communication protocols. 

Photon subtraction acting on a TMSV provides a prob- 
abilistic increase in entanglement for which a trade-off 
exists between the probability of success and the entan- 
glement of the output state. In the following, we quantify 
entanglement according to the entanglement negativity 
p9j, a monotone defined as 



J^{Pab) = 



(2) 



where ||cr|| = Tr[Vcrtcr], and p^^ is the partially trans- 
posed density matrix of the bipartite system. 
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FIG. 1. Implementation of photon subtraction: photons are 
subtracted from the input mode a by means of a high trans- 
mittivity beam splitter which reflects a small fraction of the 
input intensity towards a photon-number resolving detector 
(PNRD). The PNRD detecting an event heralds successful 
photon subtraction (30. . An approximate PNRD can be conve- 
niently implemented by multiplexed detection with avalanche 
photodiodes (APDs) ^ [32] 



In order to find the measurement operator correspond- 
ing to photon subtraction as depicted in Figure [T] we use 
the expansion of the beam splitter operator [33j: 



BaaM=e- 



-ln(c 



j)(d|ds— a'd) tan |d'ds 



(3) 



where, cos(|) = t is the transmittivity of the beam split- 
ter and sin(|) = r its reflectivity. The ancillary mode 
ttg, initially in the vacuum state, interacts with mode a 
on the beam splitter and is then projected out. Using 



by acting Mi, a on I^Ptmsv)- The norm squared of 
this state gives the probability of photon subtraction 
Ps = ('0i,o|'0i,a)- Increasing the transmittivity t of the 
beam splitter used for photon subtraction increases the 
negativity of entanglement of this state but decreases 
the probability of photon subtraction. This trade-off ap- 
pears when considering different subtraction strategies, 
as shown in Figure [2j Notably, different strategies are 
advantageous in different regimes of entanglement con- 
centration. The optimal trade-off, highlighted in gray, 
will serve as a reference to analyse the efficiency of the 
schemes considered in the following section. 
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FIG. 2. Negativity of entanglement N'{pab) as a function of 
the logarithm of the success probability ps , for different pho- 
ton subtraction strategies. These quantities are reported as 
a parametric plot when the transmittivity t of the subtract- 
ing beam splitter is varied. The curves are labeled kA/ks 
where kA^ks) is the number of photons subtracted from the 
mode A(B). For reference, the negativity of entanglement of 
some subtracted states for = 0.9 as well as a Bell state are 
marked. The optimum subtraction strategy is highlighted in 
gray. 



III. COMBINING TWO GAUSSIAN STATES 

Performing CV entanglement distillation from Gaus- 
sian states and ending with Gaussian states requires de- 
Gaussification of a large number of resources followed by 
an iterative protocol that reproduces a Gaussian state 
[27J. One Gaussification step consists of the combination 



3 



Gaussian TO 
measurement vA^ 

photon y 
subtractiot70^ 

(u dicks) ^ 



, photon 
% subtraction 
(i dicks) 



A 
B 



TMSV1 




TMSV2 



photon 
subtraction 
( V clicks) 

Gaussian 
measurement 



bg^i photon 
. subtraction 
(j clicks) 



FIG. 3. First step in the iterative protocol proposed in [27 . A 
and B share two TMSVs de-Gaussified by photon subtraction. 
The states are combined locahy on a 50:50 beam sphtter. A 
Gaussian measurement leaves two of the modes in states with 
increased entanglement. 



of pairs of identical states, employing 50:50 beam split- 
ters and vacuum projection. This linear-optical scheme 
enhances entanglement as the outputs converge towards 



a Gaussian state. While the resources required for such 
a protocol increase rapidly with the number of iterations 
the first step alone provides a feasible setup for dis- 
tillation starting with a pair of two- mode states. We 
now give analytic expressions for the four-mode states 
obtained by combining different photon-subtracted states 
on beam splitters as pictured in Figure |3] Throughout 
this section, the results of non-trace-preserving opera- 
tions are given as unnormed Hilbert space vectors whose 
norm squared equals the probability of successful imple- 
mentation. The result of subtracting i and j photons, 
respectively, from the modes of a TMSV, will be denoted 
The result of subtracting u and v from modes 
ai, a2, bi and 62, respectively, as shown in Figure [3) 
and combining the two photon-subtracted states on beam 
splitters, will be denoted |^i,j,n,v). In order to obtain 
closed form expressions, we need to presenve the symme- 
try of the states as much as possible, hence consider the 
same number of beam splitters must be inserted in the 
modes of states 1 and 2. Thus, when necessary, we con- 
sider a beam splitter inserted on one side but no photons 
detected in the reflected mode. We denote these events 
with the symbol 0. The results for various one-photon 
subtraction schemes are 



1^1,0,6,0) = ^(l^l,0)a„6, ^ l^6,o)a„6. + l^6,o)a„6, ^ l^l'O),,,^. 5 
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1^1,0,0,1) 
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1*1,0,1,0) = ^(l^2,0),„,, ® |^0,o),„fe, + I^O,o),,,fc, ® \^2fi)a.M) 

1*1,1,1,1) = Y.^XtY^^{{n + l)(n + 2) + (m + l)(m + 2)) |n, n)„^,,^ ® |m, m)„^,,^ + 



(6) 



The method employed to derive the expressions is de- 
scribed in the Appendix. We now discuss the Gaus- 
sian measurement in Figure [3] In the original proposal 
[27J, vacuum projection was used. Here, we look at a 
more general scenario, with Gaussian measurements im- 
plemented by balanced homodyne detection (i.e. a pro- 
jection on quadrature eigenstates) and 8-port homodyne 
detection (i.e. projection on coherent states). The use of 
the latter setup for vacuum projection has been discussed 
previously [28} [35j. 

The phase space distribution of a Gaussian state is 
completely described by the displacement vector (mean 
position in phase space) and the covariance matrix of 



the quadrature amplitudes. The covariance matrix of 
states conditioned on the result of a Gaussian measure- 
ments on Gaussian states does not depend on the partic- 
ular measurement result. Only the displacement of the 
obtained state is a function of the measurement result. 
Therefore, given the result and the initial state, the out- 
put can always be centered to the origin by appropriate 
post-measurement displacements ^T0\ [TT] . Therefore if 
the Gaussian measurement is performed before the sin- 
gle photon detections, for any given outcome the state 
corresponding to a null displacement vector can be ob- 
tained by classical communication and local displacement 
of modes ai, 61, a^i, as2, ^si? and 6^2, as denoted in Fig- 
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ure[3] Since the entanglement properties of states depend 
only on their covariance matrix, we can assume that the 
Gaussian projections in the setup of Figure [3] are imple- 
mented deterministically. 

For 060) 1^1,0,1,0)7 the negativity of entan- 
glement of the output is maximum if the two measured 
modes are projected on coherent states. For 1^-^ 100)' 
1^1,0,0,1) and 1^1,1,1,1), simple homodyne measurements 

- projections on eigenstates of X for a2 and of P for 62 

- are more favourable. The states obtained by measure- 
ment are proportional to the vectors |0i,j,n,v): 



k^i,o,o,o) 


= lV^l,0)ai,6i5 




^1,1,0,0) 
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h2 + (At2)2(n+l))|n,n)„^^,^ 


1^1,0,0,1) 


= ^(Atr(n + 

n 




1^1,0,1,0) 






1^1,1,1,1) 
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+ l)(n + 2) + 2+ 



+(At2)2(((At2)2 + 2)n(n + 3) + 2{\ef + 3)) |n, n) 

(7) 
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FIG. 5. Negativity of entanglement that can be obtained 
by the (1, 0, 0, 1) measurement as a function of the homodyne 
detection efficiency 77 (solid line) compared to that achieved by 
photon subtraction from a TMSV, with the same probability 
(dashed line). The transmittivity of the subtracting beam- 
splitter is chosen so that, for perfect detection, the difference 
between the values of negativity of entanglement achievable 
by the two compared strategies is maximum. 

It is then interesting to investigate the effect of lossy 
detection for this particular state. The setup is revealed 
to be fragile. Entanglement is increased only for detector 
efficiencies rj well above 0.9, as shown in Figure |5] This 
^5pould be demanding for detectors working in the pulsed 
regime. 



The probability of success and negativity of entanglement 
that can be achieved with the five measurements consid- 
ered are shown in Figure [4j In comparison with a sim- 
ple photon subtraction strategy, only the state 0,0,1) 
presents a better trade-off between increased entangle- 
ment and the success probability. 
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FIG. 4. Negativity of entanglement versus the logarithm of 
the success probability for different joint measurement strate- 
gies, yielding states given in Eq. [7| for all values of the trans- 
mittivity t of the used beam splitter. The corresponding strat- 
egy is denoted next to each line, according to the notation 
used throughout this section. The gray curve represents the 
negativity of entanglement that can be achieved by photon 
subtraction from a TMSV, as found in Figure |3] 



IV. CONCLUSIONS 

We have studied entanglement concentration by pho- 
ton subtraction from two-mode squeezed states, within 
two types of setups: (i) simple photon subtraction from 
one two-mode squeezed state and (ii) photon subtrac- 
tion from a pair of two-mode squeezed states and an 
implementation of the first step in the iterative distil- 
lation protocol, allowing a modified Gaussian measure- 
ment [27^^28]. Our results show that the second type of 
setup can, in ideal conditions, provide an advantage in 
terms of gained negativity of entanglement for a given 
success rate. However, the arrangement that allows this 
is highly sensitive to detector imperfection. Furthermore, 
a requirement for entangled states to be useful resources 
for CV quantum communication is that these states have 
a high fidelity with the closest Gaussian states [T9j. The 
best-performing two-source arrangement is not symmet- 
ric and therefore does not fall within the class of oper- 
ations that guarantee an evolution towards a Gaussian 
state. Combining our results with those in earlier work 
[36j, we are lead to conclude that meeting the require- 
ments for an improvement over photon subtraction for 
entanglement concentration, in the presence of loss, is 
likely to demand and motivate further technological de- 
velopments. 
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APPENDIX 

Here, we explain the raethod by which the expressions 
of Eq. [6] can be obtained. To this end, we give the exam- 
ple of |^o,o,o,o)- Two initial two- mode squeezed vacuum 
states are represented: 

|V^o,o)^|V^o,o) = (1-A') ^ A^+"^ |n)^(o) |n)^(o) |m)^(o) |m)^(o) 

(8) 

Reordering the modes and applying the beam-splitter 
transformation to the creation operators yields: 
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|*o,o,o,o) = Sai,a.4i,6.(l " A^) ^ A"+™ |n) |m) |n) |m) 

n,m 

^(0)n ^(0)m r(0)n r(0)m 

= (1 - a^,|:a"-b...„!L^b:.,„ |„„,B.„,^i=B. 100, 

= (1 - (^)" (^)" (^)" l^)" = 

= (1 - 5: A-+" ^.^^ E ( I' ) ( T ) ( p ) ( 7 ) (-i)'+'«r-'-'«r'i.r"-''-''?*' ioooo> = 

n,m k,l,p,q / \/ \ / 

n,m k,l,p,q ^ ^ ^ / \ / \ / 

^J(n^m- p- q)\{p + ^)! |n + m - /c - /) |/c + /) |n + m - p - g') |p + g^) ; 
Introducing iV = n + m,i^ = /c + / and P = p + g, |^o,o,o,o) = ^^/(^' ^' ^) 1^ " ^) 1^) 1^ " ^) 1^) 

„h»e = ))--'ei:e(-.)'--- (:)( ^ ^ 

(9) 



It can be shown that f{N^K^P) = 5k, p- Relabehng 
i = N — K and j = K this yields: 

1^0,0,0,0) = (l-A^)E^^^''|^)aJ^)5. \j)a. 

i,j (10) 
= |V^o,o) |^o,o) . 

This result can be confirmed by employing the Gaus- 
sian states formalism [JHS] and calculations similar to the 
one above can be performed for (non-Gaussian) photon- 
subtracted states. Functions of the type f{N^K^P) 
can be calculated, as they are simply sums of hy- 
pergeometrical terms. Let us denote the summand 
F{N^ P, n, k^p). We have (i) the summand is outside 
the summation intervals and (ii) ratios of the type: 

F{N + JAT, K + jfe, P^jQ.n^ jn, k + j/c,p + jq ' 

where the i's an the j's are integers, are rational func- 
tions. f{N^K,P) can be obtained by induction, after 
finding a recurrence relation for the sum. One way 
to get this is to find a recurrence for the summand 
F{N, P, n, /c,p), which is free of n, k and p. 
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For clarity, let us look at what this means in the case of 
a hypergeometric term in two indices: r and s (the task 
is to sum over s). For some integers / and J, one needs 
to find polynomials aij{r) that satisfy: 

I, J 

a,,,(r)P(r + j,s + z) = 0. (11) 

i,j=0 

This recurrence relation can be summed over all s. Be- 
cause the coefficients of P(r, s) in the recurrence are free 
of s and because the summation interval can be extended, 
this yields a recurrence relation for the sum, /(r): 

J 

Y a^,,(r)/(r + j) = 0. (12) 

i,j=0 

Sister Celine's algorithm [37] is a way of finding recur- 
rences free of a number of indices. The algorithm is de- 
scribed in [38j. The required induction is obtained com- 
putationally, using the code MultiSum [39j, written in 
Mathematica [40]. 



